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The Configuration Based Aero dynamics (CBAERO) software package is used 
to predict the convective and radiative heating environments for the Crew 
Exploration Vehicle (CEV). A limited number of high fidelity CFD solutions are 
used to ‘anchor’ the engineering level estimates obtained using CBAERO. 

I. Introduction 

The Configuration Based Aero dynamics (CBAERO) software package, described in reference 1, is an 
engineering level aero-thermodynamics tool for predicting the aerodynamic and aero -thermodynamic environments 
of general vehicle configurations. CBAERO is based on accepted engineering level analysis for the estimation of 
both aerodynamic and aerothermal environments. In this paper a new anchoring process is described wherein high 
fidelity CFD results are used to ‘anchor’ the engineering level analyses. The anchoring process permits a few, select, 
CFD solutions to be used to beyond the specific flight condition they were original run at, allowing for the cost 
effective use of high fidelity, and computationally expensive, CFD solutions early in the design process when the 
vehicle trajectories are often in a constant state of change. 

II. CBAERO 

The Configuration Based Aero dynamics (CBAERO) software package, described in reference 1, is an engineering 
level aero-thermodynamics tool for predicting the aerodynamic and aero-thermodynamic environments of general 
vehicle configurations. For the present work CBAERO is used as the basis of the engineering level aerothermal 
analysis upon which corrections, based on CFD, are applied. 

A. Geometry Definition 

CBAERO makes use of an unstructured surface grid of triangles to define the Outer Mold Line (OML) of the 
vehicle configuration. No volume mesh is required. Figure 1 presents the unstructured surface mesh used to define 
the Crew Exploration Vehicle (CEV) geometry. The surface mesh for the Vi body contains 5,396 triangles. 

B. Subsonic Mach Numbers 

For the subsonic Mach number range CBAERO makes use of a fast, multi-pole, unstructured panel code 
formulation. The linear, integral formulation for subsonic, incompressible flow is solved using a fast multi-pole 
method, sometimes generically referred to as a fast tree method [3, 4]. An octree data structure is used to define the 
spatial distribution of the surface triangles (panels), and their relative proximity to one another. Using the octree 
information, approximate, multi-pole, expansions are built for entire regions of the of the solution domain. The 
solution of the linear system of equations is then performed in ~N*Log N time, [3, 4] using a preconditioned 
GMRES algorithm. 

C. Super and Hypersonic Mach Numbers 

For super through hypersonic Mach numbers the inviscid solution is based on independent panel methods, such 
as Modified-Newtonian, Tangent Cone, or Tangent Wedge formulations. All results presented here made use of the 
Modified-Newtonian formulation. Figure 2 depicts the surface pressure contours, using a Modified-Newtonian 
formulation, on the CEV configuration. 
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The independent panel methods provide the pressure at each surface triangle. The entropy on each triangle is 
taken equal to the post normal shock entropy. The known value of pressure and entropy at each panel fixes the 
thermodynamic state at the edge of the boundary layer. A table lookup procedure with data created using the Gordon 
and McBride (CEA) code, of reference [8] is used to calculate the remaining thermodynamic properties such as 
temperature and enthalpy. 



Fig. 1 Unstructured Triangulated Surface for the CEV, containing 5396 triangles. 



Fig. 2 Mach 32.2, 157.0° Angle of Attack, CBAERO Surface Pressure Field. 
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D. Surface Streamline Pattern 


For subsonic flows the panel solution provides the local surface velocity. For supersonic and hypersonic flows 
the following approximation for the surface velocity field is used: 


surface 


-hxV^xh 


Here h is the local surface normal, and V CX) is the free stream velocity vector. This approximation is used to define 
the entire surface velocity flow field. 

A robust integration algorithm is used to trace the streamline pattern over the entire vehicle. The algorithm 
calculates the streamline pattern in reverse, starting at the rear of the vehicle and integrating forward, until the 
stagnation point is reached. The algorithm continues to generate streamlines until every node/triangle in the mesh has 
a nearby streamline from which the distance to the stagnation point or an attachment line can be calculated. 

The method of reference [2] is used to locate the attachment lines. The method makes use of topology 
considerations to locate attachment and separation lines. Here, only the attachments lines are of interest. The running 
lengths from either the stagnation point or the attachment lines are then calculated. These running lengths are used in 
the acreage heating models in the calculation of the local Reynolds number. Figure 3 depicts the streamlines on the 
CEV, calculated using the above procedure for a free stream Mach number of 32.2 and an angle of attack of 157°. 


E. Blanked or Shadowed Regions 


CBAERO automatically calculates those regions that are ‘blanked’ or ‘shadowed’ by other windward facing 
surfaces. A sorting algorithm is used to quickly locate any windward facing surfaces that are ‘blanked’ by other 
windward facing surfaces. The pressure coefficient on blanked surfaces is set to an empirically derived value which 
tends to zero at high Mach numbers. 



Fig. 3 Mach 32.2, 157.0° Angle of Attack, CBAERO Surface Streamline Pattern 


F. Stagnation Point Heating 

The convective heating correlation of Tauber [5] is used to predict the convective heating at the stagnation 
point. The correlation is driven by the stagnation edge conditions, as well as the local radius. The radius is estimated 
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by fitting a general bi-cubic surface at the location of the stagnation point. Sufficient neighboring triangles are used 
to calculate a best fit surface, in a least squares sense. The primary curvatures are then calculated directly from this 
cubic surface. An average of the principle curvatures is used to determine the local radius. For a perfectly spherical 
stagnation point both principle curvatures will be identical and the local radius estimate will reduce to the expected 
radius of the sphere. 

G. Attachment Line Heating 

As with the stagnation point, the convective heating correlation of Tauber [5] is used to predict the convective 
heating along the attachment lines. The correlation is driven by the stagnation edge conditions, as well as the local 
radius. The radius is estimated by fitting a general bi-cubic surface at the location of each of the previously identified 
triangles lying on the attachment line. Sufficient neighboring triangles are used to calculate a best fit surface, in a 
least squares sense. 

When fitting the bi-cubic patch, the two local surface coordinate directions are aligned parallel and normal to 
the local surface velocity. The third coordinate direction remains normal to the surface. The primary curvatures are 
then calculated directly from this cubic surface. The surface curvature in the direction normal to the local surface 
velocity is used to determine the local radius. For a simple swept cylinder this formulation would return the expected 
radius of the cylinder. This radius, along with an appropriate sweep correction is used in Tauber’s correlation to 
predict the attachment line convective heating. 

H. Radiative Heating 

The Tauber-Sutton radiative heating correlation of reference [9] is used to predict the stagnation and off 
stagnation radiative heating environments. The correlation is of the form 

Cp A R‘-f(V) 

Here p is the free stream density, R is the local radius at the stagnation point, V is the velocity normal to the local 
shock surface, f(V) is the intensity function, and A, B, and C are constants. The local radius, as determined for the 
convective heating model, is again used. The off stagnation point distribution is scaled from the stagnation point 
value with an empirically derived function that is a function of the local impact angle. 

I. Acreage Heating 

The acreage heating model is based on either laminar or turbulent flat plate reference enthalpy methods for the 
skin friction, and Reynold’s analogy to estimate the Stanton number and hence the local heating rate. A simple 
Mangier correction is used to adjust the two-dimensional estimates. The running lengths, calculated using the 
streamline patterns, are used to estimate the local Reynolds number. For both laminar and turbulent flows a reference 
temperature, based on Eckert’s method, is calculated: 

T ref = 0.5 • (r e + rj + o.i i • • O' - 1) • m] ■ T e 

The recovery factor, R f , is calculated as the square root of the Prandtl number for laminar flows, and the cube 
root of the Prandtl number for turbulent flows. The Prandtl number is evaluated at the reference temperature. For 
laminar flows the skin friction coefficient and Stanton number are calculated as: 

C f = 0.664 / a /R d~M~ f ■ T e / T ref 
S t - 0.5 • C f ■ Pr“ 2/3 

For turbulent flows the relations for skin friction and Stanton number are: 

C f - 0-37 /[log 10 (Re/ M f )f 584 . T e / T ref 
S t = 0.5 • C f • Pr -2/3 

In the above relations, M f is the Mangier factor. For laminar flows, M f is equal to 3 and for turbulent flows it is 
set to 2. The enthalpy based film coefficient, recovery enthalpy, and convective heating are defined as: 

C H ~ Pref * U e * St 
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H rec ~H e + 0.5 • R f ■ U 2 e 

Qconv = C H •( H rec ~ H J 

In the above relations, H e is the edge static enthalpy, H rec is the recovery enthalpy, U e the edge velocity, and H w 
is the wall enthalpy. 

The Reynolds number used above is based on the geometric running lengths estimated using the Newtonian 
surface velocity approximation for super or hypersonic flows or the fast-panel solution for subsonic flows. On 
windward surfaces the running length is taken directly as the geometric value, on leeward surface a weighted running 
length is used. The process is similar to that used in reference [17]. For both laminar and turbulent flows, the 
weighted reference length, S w , is calculated as 

p eJ u e U e S w = J p eJ u e U e dS 

0 

The integration is along a streamline. A smooth transition, from the geometric to the weighted running length is used 
as the streamline passes from the windward to leeward surface. Finally, it is noted that in the calculation of the local 
Reynolds number the reference density and viscosity are used. 

J. Transition 

CBAERO can be run fully laminar, fully turbulent, or transitional. For transitional flow a simple model, based 
on Re^/M edge is used to determine the onset of transition. The transition length is set equal to the laminar run 
length, with a linear transition from laminar to fully turbulent flow. 

K. Free Molecular and Rarefied Flow 

CBAERO does include both a free molecular flow model [18] and bridging function methodology [19] for 
rarefied flows. However, for the present work only the above continuum models were used for the databases created 
for the CEV thus far. 

L. Catalytic Surface Model 

To date, all the aerothermal databases for the CEV have been run fully catalytic. However, CBAERO does 
include partially catalytic heating models. The partially catalytic heating models developed for CBAERO incorporate 
mass transfer coefficients and global surface reaction rates for a chemically frozen boundary layer [14]. The wall 
surface reaction rates are determined by the speeds at which oxygen and nitrogen atoms move towards the wall, 
which in turn depend on the wall temperatures and the atomic masses of Oxygen and Nitrogen. The wall reaction rate 
constants are calculated separately for Oxygen and Nitrogen. The recombination coefficients are normalized wall 
reaction rates which are found experimentally to be strong functions of surface temperatures and material 
characteristics [7]. 

A similarity variable, (3 , represents the velocity gradient at the stagnation point for a hemisphere. Calculated 
from Newtonian theory, it depends on the free stream velocity, nose radius and the free stream and boundary layer 
edge densities, 

P = ( V J R m J-p- 0 *P-/P^e 

For the acreage heating points it can be demonstrated that the similarity variables are the velocity at the 
boundary layer edge and the running length. The similarity variable is taken to be of the form 

j3 = C-{VJX)-f{S/R nose ) 

Where C is a constant dependent on the form of the laminar flat plate heating model, X is the local running length, 
and f(S,R) is an exponentially decaying function fitted to smoothly match the flat plate catalytic surface model with 
the stagnation catalytic surface model. 

The fraction of atoms that reach the wall and recombine is represented by (j) 0 N which is calculated for both 
Oxygen and Nitrogen. The convective heat transfer rate for a fully catalytic surface can be written as: 
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a FC 

“ convective 


= H r (I F R c -I wall ) 


Here ‘FC’ stands for ‘Fully Catalytic’. 

FC 

It should be noted that the recovery enthalpy in a fully catalytic formulation, I R , contains the entire energy of 
the O and N atoms which have dissociated at the boundary layer edge. For the partially catalytic case, however, a 
fraction of these atoms, (l.O — ^ 0N ), do not recombine on the surface. Therefore, the fully catalytic recovery 

PC 

enthalpy is replaced by the partially catalytic recovery enthalpy, / R , by subtracting the enthalpies of the atoms 
which do not recombine: 


C c = 

0,N 

Here ‘PC’ stands for ‘Partially Catalytic’, represents the mass fractions of Oxygen and Nitrogen atoms at the 

boundary layer edge, and /q N are their heats of dissociation or recombination. These mass fractions are determined 

by the chemical equilibrium program CEA [8] at the boundary layer edge conditions and stored in a tabular fashion 
for fast lookup within CBAERO. 


M. Surface Energy Balance 


For any triangle on the surface an energy balance for a non-conducting wall is formulated. In general, the energy 
radiated from the surface must equal the sum of the convective and incident shock-radiative heating: 


+ <1 radiative ~ G ' 8 ' C 


wall 


Since the convective heating depends on the final wall temperature this non-linear relationship must be solved 
iteratively at each triangle. 


III. CFD 

To judge the accuracy of CBAERO, comparisons have been made to the NASA Ames CFD codes DPLR and 
NEQAIR. DPLR [10] is a parallel multi-block finite volume code that solves the reacting Navier-Stokes equations 
including finite-rate chemistry and the effects of thermal non-equilibrium. Based on a finite-volume implementation 
of a modified Steger- Warming flux splitting method with Gauss-Seidel line relaxation, DPLR is a truly parallel, 
multi-block codeNominal third-order spatial accuracy, for inviscid fluxes, is obtained via MUSCL extrapolation 
coupled with a minmod limiter. Viscous fluxes are computed to second-order spatial accuracy using central 
differencing. 

NEQAIR [11] is a line-by-line spectroscopy and one-dimensional radiation transport tool. The code computes 
radiation transport (from the free stream to the vehicle surface) along straight lines of sight using a tangent-slab 
approximation. The data at points on a line of sight are obtained through interpolation from a volume solution 
computed using DPLR, i.e., radiative heating is obtained a posteriori using converged flow solutions, essentially 
assuming the flow and radiative processes are uncoupled. NEQAIR constructs accurate spectra at every point on a 
line of sight, computes the integrated radiative intensity from the spectral distribution, and transports the energy to 
the next point on the line of sight through one-dimensional transport. 

A. Comparison of CBAERO to CFD 

In figures 4 through 8 CBAERO results for a Mach 28.4, 152° angle of attack, laminar case are compared with 
results from the NASA Ames CFD codes DPLR and NEQAIR [10, 11]. For this case, the CBAERO results took 3.7 
seconds of CPU time on a 2.26 Pentium M laptop. 

The CBAERO pressure distribution presented in figure 4 compares favorably with the DPLR result. Similarly 
the CBAERO prediction for convective heating is presented in figure 5. The CBAERO results compare nicely with 
the DPLR results, slightly under predicting the peak shoulder heating by approximately 10%. 

Figure 6 compares the CBAERO temperature distribution with the DPLR highlighting the details of the leeward 
heating environment. Once again, CBAERO does a reasonable job of predicting the leeward heating environment. 
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Figure 7 compares the CBAERO and DPLR predictions for the laminar shear magnitude. CBAERO clearly 
under predicts the peak shear stress near the shoulder, but otherwise predicts the general shape of the shear stress 
distribution. 

Figure 8 compares the CBAERO and NEQAIR [11] results for the radiative heating. The peak radiative heating 
magnitude compares favorably with the NEQAIR results, with the peak slightly offset from the NEQAIR results. The 
over all shape of the radiative heating distribution is generally captured. 

Figures 9 and 10 compare the surface distribution of the CBAERO laminar convective heating and radiative 
heating predictions with DPLR. As was evident in the centerline plots, CBAERO tends to under predict the laminar 
convective heating near the shoulder and over predict across the heat shield. For the radiative heating, CBAERO 
captures the correct maximum value but places the peak closer to the shoulder than does NEQAIR. Over the heat 
shield, CBAERO tends to slightly under predict the radiative heating. 


Centerline Slice for M: 28.409700. Rho: 0.000350. T: 249.770004. AoA: 152.000000 Pressure vs. Z 

CFD Case: 29 

CFO File: . . /CF D_D ata/dplrji i_speed/cev .basel ine.9kms.a28 . la m 



Z. m 


Fig. 4 Mach 28.4, 152.0° Angle of Attack, Centerline Pressure Distribution 
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Fig. 5 Mach 28.4, 152.0° Angle of Attack, Centerline Laminar Convective Heating Distribution 


Centerline Slice for M: 28.409700, Rho: Q.C0035CL T: 249 7700-34, AoA: 152.000000 Temperature vs. Z 

CFD Case: 29 

C— C Fils: ...'CFD_Data/dp r_h_speed/cev.base ire.9krns.a28 lam 



Z. m 

Fig. 6 Mach 28.4, 152.0° Angle of Attack, Centerline Laminar Temperature Distribution 
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Fig. 7 Mach 28.4, 152.0° Angle of Attack, Centerline Laminar Shear Stress Magnitude Distribution 


Centerline Slice for M: 28.409700, Rho: 0.000350, T: 249.770004, AoA: 152.000000 Radiative Qdot vs. Z 

CFO Case: 29 

CFO File: ../CFD_Data/dplr_hi_speed/cev.baselin€.9kms.a28.lam 
CFD File: . ./CFD_Data/dpii_hi_speed/neqair_resuits/T Cl Siam. Rad 



-3 -2 -1 0 12 3 


z, m 

Fig. 8 Mach 28.4, 152.0° Angle of Attack, Centerline Radiative Heating Distribution 
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Fig. 9 Mach 28.4, 152.0° Angle of Attack, Surface Laminar Convective Heating - CBAERO solution compared with DPLR 
results. 



Fig. 10 Mach 28.4, 152.0° Angle of Attack, Surface Radiative Heating - CBAERO solution compared with 
DPLR results. 
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IV. The Anchoring Process 


The CPU requirements for CBAERO are minimal, typically requiring less than 10 seconds per solution on the 
current configuration and mesh. This is significantly less than that required for a single CFD solution from either 
DPLR [10] or LAURA [15, 16]. To maximize the return on the computational resources invested in the DPLR and 
LAURA computations, the CLD solutions are used to anchor the CBAERO results. 

At the limited number of free stream conditions (Mach, Dynamic Pressure, and Angle of Attack) at which CFD 
solutions are available the differences in the CFD and CBAERO solutions are compared. A companion piece of 
software is used to generate correction factors for the following CBAERO quantities: 1) pressure, 2) laminar 
convective heating, 3) turbulent convective heating, 4) laminar skin friction magnitude, 5) turbulent skin friction 
magnitude, 6) radiative heating, and 7) the boundary layer edge velocity components. The correction factors are 
calculated for each surface triangle for each CFD Mach-Q- Alpha condition. The corrections factors are then linearly 
interpolated in Mach-Q-Alpha, Mach-Q, Mach-Alpha, or Q-Alpha space as defined by the user. The linearly 
interpolated corrections are then used to anchor CBAERO solutions at free stream conditions for which no CFD 
solution is available. 

A. Pressure 

For the pressure corrections, the assumption is made that the laminar and turbulent pressure distributions are 
similar enough that only an anchoring using the laminar CFD solution is required. The pressure ratio correction is of 
the form 


Pratio — ‘ 


CFD Pressure 
CBAERO Pressure 


Here the pressures are the local pressures for each CBAERO surface triangle, obtained by interpolation from the 
original CFD surface mesh. The value of P ratio is limited to values greater than 0 and less than 20. 

B. Convective Heating 

For the laminar and convective heating the correction ratios are of the form 

CFD Laminar Q 


Clai 


CBAERO Laminar Q 


Ctrb . = - 


CFD Turbulent Q 
CBAERO Turbulent Q 


Here Clam ratio and Ctrb ratio refer to the correction ratios applied to the CBAERO estimates for the laminar and 
turbulent heat transfer coefficients. The values of the laminar and turbulent heating correction ratios are limited to 
values greater than zero and less than 10. Solving for Clam ratio and Ctrb ratio is an iterative process due to the nonlinear 
nature of the surface energy balance equation. 

C. Shear Forces 

The anchoring of the shear force magnitude requires a minor modification to the ratio formulation. Since shear 
force may go to zero, the correction ratio takes the form 


Cflai 


CFD Shear Magnitude 
CBAERO Shear Magnitude + G 
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Cftrb 


ratio 


CFD Shear Magnitude 
CB AERO Shear Magnitude + G 


Here S is a small number, taken to be 0.25 Pascals. The values of Cfiam or Cftrb are limited to values greater 

7 ratio ratio 0 

than 0 and less than 10. 


D. Radiative Heating 

The anchoring for radiative heating is similar to that for the shear magnitude 


Q® ratio = 


CFD Radiative Q 
CB AERO Radiative Q + £ 


Once again, £ is a small number, taken to be 0.25 Watts/m A 2. The value of Qs ratio is limited to values greater than 
0 and less than 10. 

E. Surface Streamlines 

The correction ratios calculated above for the laminar and turbulent skin friction account only for differences in 
the magnitude of the CBAERO and DPLR results. To further improve the CBAERO model, the differences in the 
CBAERO predicted surface streamline pattern, and those calculated by DPLR are accounted for. The CBAERO 
approximation for the local surface (or edge) velocity is simply 

^surface = nXV^Xh 

Based on the x, y, and z components of the shear forces calculated by DPLR a correction factor can be calculated 
such that the adjusted surface (edge) velocity used by CBAERO accurately recovers the DPLR streamlines (based on 
shear). Unlike the previous correction ratios, the surface velocity corrections are deltas. 

F. Interpolating the Corrections to New Flight Conditions 

The correction ratios for pressure, convective heating, shear force magnitude, and radiative heating are 
calculated at each available CFD solution. The CFD anchoring ‘space’ is defined by a small number of CFD 
solutions in Mach - Dynamic Pressure - Alpha space. To apply these corrective ratios at any Mach-Q- Alpha 
condition an interpolation process must be defined. Currently, the CFD points are meshed in two or three dimensions 
using an unstructured tessellation of triangles or tetrahedra. The two or three dimensional unstructured mesh is 
created to satisfy the Delaunay criterion and thus enclose the convex hull of the CFD points. Figure 11 shows a 
hypothetical mesh in Mach-Q- Alpha space defined by a single tetrahedron. A similarly hypothetical and simple mesh 
in two dimensions would contain three nodes and a single triangle. 

Linear interpolation is performed on the unstructured mesh to calculate the correction ratios for Mach-Q- Alpha 
condition that lies within the convex hull of the CFD points. For a Mach-Q- Alpha condition outside the convex hull 
value returned from the interpolation (extrapolation) is ‘clipped’ to the edge of the ‘table’. Hence extrapolation of 
the corrections is explicitly forbidden. 
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D. Testing of the Process 

As a test, of the anchoring procedure, approximately V 2 of the CFD solutions available for anchoring are 
removed from the ‘anchoring list’. The remaining CFD solutions are then used to create the anchoring database as 
described above. CBAERO is then run, using the anchoring corrections, at the reserved CFD conditions and 
comparisons are made. The total list of available CFD solutions is presented in Table 1. The first 14 CFD solutions 
in Table 1, highlighted in red, are used in the anchoring procedure. These first 14 CFD solutions are those originally 
generated during NASA’s Exploration Systems Architecture Study (ESAS). Figures 12 and 13 present the Mach- 
Alpha tessellation of the CFD point set used for this anchoring test. 

Figure 12 shows the two dimensional tessellation used to interpolate the correction ratios in Mach- Alpha space. 
Figure 13 shows the true distribution of the CFD anchor points in the complete three dimensional Mach-Alpha-Q 
space in which they reside. The uneven and sparse nature of the dynamic pressure variance (Q) of the CFD anchor 
points was the driver for restricting the interpolation to two dimensions (Mach and Alpha). 

Figures 14 through 18 present the anchored CBAERO results, using the first 14 CFD solutions, compared to 
DPLR for the Mach 28.4, 152° case presented earlier. These results should be compared with figures 4 through 8. 

The anchored CBAERO pressure distribution along the centerline, shown in figure 14, matches the DPLR 
results very nicely. Some differences are noticeable on the leeward side, where the anchoring ratios were likely 
clipped. Figure 15 compares the anchored CBAERO results for the laminar convective heating with DPLR. Once 
again, the anchored results compare very nicely with the exact DPLR results, even on the leeward side. Figure 16 
compares the anchored CBAERO temperature distribution, highlighting the leeward heating environment. Once 
again, the anchored CBAERO results compare very favorably with the exact DPLR solution. Figure 17 presents the 
anchored CBAERO laminar shear magnitude results. The anchored results miss the peak shear magnitude on the 
shoulder by approximately 20%, however this is an improvement over the unanchored results which missed the shear 
magnitude by a factor of 2.5. Finally, figure 18 compares the anchored CBAERO radiative heating prediction with 
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the NEQAIR results. The anchored CBAERO more closely matches the NEQAIR radiative heating result in both 
magnitude and general shape. 

Figures 19 and 20 compare the surface distribution on the anchored CBAERO laminar convective and radiative 
heating with the DPLR and NEQAIR results. These anchored CBAERO results can be compared with those in 
figures 9 and 10. 

V. Final Aerothermal Database Creation 

The final anchored CBAERO laminar and turbulent databases use all of the available CFD anchor points from 
DPLR, LAURA, and NEQAIR. Figure 21 shows the two dimensional anchoring mesh used in Mach- Alpha space. 
Figure 22 shows the true distribution of the anchoring points in the complete Mach-Alpha-Q space. Once again, due 
to the uneven and sparse nature of the CFD anchor points in the Q-dimension the interpolation of the anchoring 
ratios was limited to Mach- Alpha space. 

Table 2 lists the Mach, Q, and Alpha conditions covered in the anchored laminar and turbulent databases. The 
anchored databases include 25 Mach numbers, 19 dynamic pressures, and 11 angles of attack for a total of 5225 
cases. Figures 23 and 24 show the distribution of the database points with the CFD anchoring points overlaid. Both 
databases were run on NASA’s Columbia super computer using 55 processors. The anchored laminar database 
required 10 wall clock minutes, while the anchored turbulent database required 15 wall clock minutes. These timings 
average to approximately 6 to 10 seconds per solution. 

Figures 25 through 31 compare the anchored CBAERO solution with the DPLR and NEQAIR results. Since this 
Mach 28.4, 152° case is now part of the anchoring CFD cases the results shown in these figures should, within the 
limitations imposed by the ratio limiters, recover the CFD solution. This is true for all the quantities for which 
anchoring is applied. The largest discrepancies appear in the maximum shear force near the shoulder, and near the 
point of minimum radiative heating on the leeward shoulder. 

VI. Conclusions and Future Work 

The CBAERO software package has been used to predict the convective and radiative heating environments for 
the Crew Exploration Vehicle (CEV). A limited number of high fidelity CFD solutions are used to ‘anchor’ the 
engineering level estimates of pressure, convective heating, radiative heating, and shear stress. This mixture of fast 
engineering level analysis, anchored by high fidelity CFD, provides for a very flexible tool during conceptual and 
preliminary design. The used of CBAERO and the anchoring process increases the value of the available CFD 
solutions beyond the single flight condition at which they were original run. The anchoring processed described here 
is currently being used to support the CEV program, and on going work is aimed at extending the anchoring process 
to include additional boundary layer edge information extracted from the CFD solutions, as well as extending the 
process to include rarefied flow solutions from NASA’s Direct Simulation Monte Carlo (DSMC) code DAC [20]. 

References 

[1] Kinney, D.J., “Aero -Thermodynamics for Conceptual Design”, AIAA Paper 13382, January 2004, Reno 
Nevada 

[2] Ken wright, David N., “Feature Extraction of Separation and Attachment Lines”, IEEE Transactions on 
Visualization and Computer Graphics, Vol. 5, No. 2, April-June 1999 

[3] Barnes, J.E., “N-Body Models of Collisionless Systems”, volume 1, chapter 8, Springer- Verlag, 1996 

[4] Barnes, J.E., and Hut, P., “A hierarchical 0(N Log N) force calculation algorithm”, Nature, 324(4):446-449, 
December 1986 

[5] Tauber, Michael E., “A Review of High-Speed Convective, Heat-Transfer Computation Methods”, NASA TP 
2914, 1989 

[6] Goulard, R., “Recombination Rates in Hypersonic Stagnation Heat Transfer”, in An Introduction to 
Hypersonic Aerodynamics , presented at the ARS Annual Meeting, NY, December, 1957. 

[7] Stewart, David A., "Surface Catalysis and Characterization of Proposed Candidate TPS for Access-to-Space 
Vehicles", NASA Technical Memorandum, 1 12206, July 1997. 

[8] McBride, B. J. and Gordon, S., "Computer Program for Calculation of Complex Chemical Equilibrium 
Compositions and Applications II. User's Manual and Program Description", NASA RP-1311-P2, June 1996 

[9] Tauber, M.E. and Sutton, K., “Stagnation-Point Radiative Heating Relations for Earth and Mars Entries”, 
Journal of Spacecraft Volume 28, Number 1, January-February 1991, Pages 40-42 


14 



[10] Wright, M.J., Prabhu, D.K., and Martinez, E.R., “Analysis of Afterbody Heating Rates on the Apollo 
Command Modules, Parti: AS-202”, AIAA Paper 2004-2456, 37th AIAA Thermophysics Conferenence 

[11] Whiting, E.E., Yen, L., Arnold, J.O., and Paterson, J.A, “NEQAIR96”, Nonequilibrium and Equilibrium 
Radiative Transport and Spectra Program: User’s Manual”, NASA RP-1389, Dec. 1996 

[12] Lee, D.B. and Goodrich, W.D., “Heat Transfer Rate and Prssure Measurements Obtained during Apollo 
Orbital Entries”, NASA TN D=6028, Oct. 1970 

[13] Cornette, E.S., “Forebody Temperatures and Calorimeter Rates Measured During Project Fire II Reentry at 
11.35 km/s”, NASA TM X-1305, Nov. 1966 

[14] Adelman, H. “HAVOC Aerothermal Heating (Partially Catalytic Surfaces)”. ELORET Task 8, Contract No. 
NAS2-0062, Report No 2, February 2001. 

[15] Gnoffo, P. A., Gupta, R. N., and Shinn, J., “Conservation Equations and Physical Models for Hypersonic Air 
Flows in Thermal and Chemical Nonequilibrium,” NASA TP 2867, 1989 

[16] Gnoffo, P. A., “An Upwind-Biased, Point-Implicit Relaxation Algorithm for Viscous, Compressible Perfect- 
Gas Flows,” NASA TP 2953, 1990 

[17] Zoby, E.V, Moss, J.N., and Sutton, K. “Approximate Convective-Heating Equations for Hypersonic Flows”, 
Journal of Spacecraft and Rockets, Vol. 18: pp. 64-70, January-February 1981. 

[18] Bird, G.A., Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Clarendon Press, Oxford 1994. 

[19] Wilmoth, R.G. and Blanchard, R.C. “Rarefied Transitional Bridging of Blunt Body Aerodynamics”; 21st 
International Symposium on Rarefied Gas Dynamics 

[20] LeBeau, G.J. “A Parallel Implementation of the Direct Simulation Monte Carlo Method,” Computer Methods 
in Applied Mechanics and Engineering Parallel Computational Methods for Flow Simulation and Modeling, 
Vol 174, pp 319-337. AIAA Paper 2006-3577, June 2006. 


15 



Mach 

Density 

Temp 

Alpha 

Lam 

Trb 

Rad 


kg/m A 3 

°K 

deg 

y/n 

y/n 

y/n 

ESAS Study Runs 

6 

0.0040000 

250.00 

147.00 

Yes 

Yes 

No 

6 

0.0040000 

250.00 

157.00 

Yes 

Yes 

No 

10 

0.0020000 

250.00 

147.00 

Yes 

Yes 

No 

10 

0.0020000 

250.00 

157.00 

Yes 

Yes 

No 

18 

0.0004000 

250.00 

147.00 

Yes 

Yes 

No 

18 

0.0004000 

250.00 

157.00 

Yes 

Yes 

No 

21.44 

0.0016271 

267.97 

153.00 

Yes 

Yes 

No 

21.44 

0.0016271 

267.97 

152.50 

Yes 

Yes 

No 

25 

0.0000400 

250.00 

147.00 

Yes 

Yes 

No 

25 

0.0000400 

250.00 

157.00 

Yes 

Yes 

No 

28.93 

0.0006269 

263.08 

153.00 

Yes 

Yes 

Yes 

29.54 

0.0005378 

259.51 

153.00 

Yes 

Yes 

Yes 

32.2 

0.0003396 

248.89 

147.00 

Yes 

Yes 

Yes 

32.2 

0.0003396 

248.89 

157.00 

Yes 

Yes 

Yes 

DPLR LEO Runs 






23.84 

0.0001632 

233.29 

150.00 

Yes 

Yes 

No 

23.84 

0.0001632 

233.29 

152.00 

Yes 

Yes 

No 

23.84 

0.0001632 

233.29 

154.00 

Yes 

Yes 

No 

23.621 

0.0001412 

230.30 

150.00 

Yes 

Yes 

No 

23.621 

0.0001412 

230.30 

154.00 

Yes 

Yes 

No 

16.78 

0.0004073 

253.15 

150.00 

Yes 

Yes 

No 

16.78 

0.0004073 

253.15 

152.00 

Yes 

Yes 

No 

16.78 

0.0004073 

253.15 

154.00 

Yes 

Yes 

No 

9.63 

0.0064270 

241 .50 

150.00 

Yes 

Yes 

No 

23.82 

0.0001412 

230.30 

152.84 

Yes 

Yes 

No 

LAURA LEO Runs 

9.26 

0.0022590 

261 .40 

150.00 

Yes 

Yes 

No 

9.26 

0.0022590 

261 .40 

152.00 

Yes 

Yes 

No 

9.26 

0.0022590 

261 .40 

154.00 

Yes 

Yes 

No 

20.4 

0.0005045 

257.98 

152.00 

Yes 

Yes 

No 

23.82 

0.0001412 

230.30 

152.84 

Yes 

Yes 

No 

DPLR High Speed Cases 

28.4097 

0.0003505 

249.77 

152.00 

Yes 

Yes 

Yes 

35.78605 

0.0000624 

214.26 

152.00 

Yes 

Yes 

Yes 


Table. 1 CFD anchoring points at which DPLR, LAURA, and/or NEQAIR solutions were available. The points in red 
indicate those CFD solutions used in the testing of the anchoring process. The LAURA solution in blue is a ‘duplicate’ of a 
DPLR solution. All CFD solutions are based on a 5.5m diameter CEV. 
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Fig. 13 
shown. 
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Centerline Slice for M: 28.409700, Rho: 0.000350, T: 249.770004. AoA: 152.000000 Pressure vs. Z 

CFD Case: 29 

CF D File: . iCF D_Data/dpl r_h i_speed/cev. baseline .9kms.a2B. lam 



-3 -2 -1 0 1 2 3 

Z, m 

Fig. 14 Mach 28.4, 152.0° Angle of Attack, Centerline Pressure Distribution - Anchored CBAERO solution compared with 
DPLR results. Only Vi of the available CFD points were used in the anchoring process, leaving the remaining CFD points, 
such as this one, for use as test cases. 


CenterlineS ce for M 28.409700 Rho 0 C0C35C T 249.770004 AoA 1 52. 0CHKMK) Convective Qdot vs Z 

CFD Case: 29 

C F D File : . . /CF D_Data/d pi r_h i_speed/cev . base! i ne .9kms„a28 . la rn 



-3 -2 -1 0 12 3 

z__ 

Fig. 15 Mach 28.4, 152.0° Angle of Attack, Centerline Laminar Convective Heating Distribution - Anchored CBAERO 
solution compared with DPLR results. Only Vi of the available CFD points were used in the anchoring process, leaving the 
remaining CFD points, such as this one, for use as test cases. 
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Centerline SI ce for M: 28.409700, Rho: 0.CCC35C. T 249 773034. AoA: 152.0 DOC Q0 Temperature vs. Z 

CFD Case: 29 

CFD File /CFD_Data/dp! ir__hijspeed/cev base ire 9kms a28 am 



Z, m 

Fig. 16 Mach 28.4, 152.0° Angle of Attack, Laminar Temperature Distribution - Anchored CBAERO solution compared 
with DPLR results. Only Vi of the available CFD points were used in the anchoring process, leaving the remaining CFD 
points, such as this one, for use as test cases. 


Center! ne Slice for 1/ 28 .4C97CC. Rho: 0.000350. 1:249.770004 AoA 152.333333 Shear Magntiuce vs Z 

CFD Case: 29 

CFD Rile: ../CFD_Data/dplr_hi_speeclfcev.baseiirie.9kiTis.a28.lam 



Z, rn 

Fig. 17 Mach 28.4, 152.0° Angle of Attack, Centerline Laminar Shear Stress Magnitude - Anchored CBAERO solution 
compared with DPLR results. Only Vi of the available CFD points were used in the anchoring process, leaving the 
remaining CFD points, such as this one, for use as test cases. 
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Centerline Slice for M: 28.409700, Rho: 0.000350, T: 249.770004, AoA: 152.000000 Radiative Qdot vs. Z 

CFD Case: 29 

CFD File: . ./CFD_Data/dpIr_hi_speed/cev.baselme.9kms.a28.lam 
CFD File: /C F D_D ata/d p lr_h !_s peed/n eq air_resu Its/TC 1 5 La m . Rad 



z. m 

Fig. 18 Mach 28.4, 152.0° Angle of Attack, Centerline Radiative Heating - Anchored CBAERO solution compared with 
NEQAIR results. Only Vi of the available CFD points were used in the anchoring process, leaving the remaining CFD 
points, such as this one, for use as test cases. 


Convective Qdojv^w/CM A 2 
307.00000 


230.2625] 
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76.78750 
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0.05000 



CBAER 


DPL 


Fig. 19 Mach 28.4, 152.0° Angle of Attack, Surface Laminar Convective Heating - Anchored CBAERO solution compared 
with DPLR results. Only Vi of the available CFD points were used in the anchoring process, leaving the remaining CFD 
points, such as this one, for use as test cases 
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Fig. 20 Mach 28.4, 152.0° Angle of Attack, Surface Radiative Heating - Anchored CBAERO solution compared with 
DPLR results. Only Vi of the available CFD points were used in the anchoring process, leaving the remaining CFD points, 
such as this one, for use as test cases. 



Fig. 21 Two dimensional tessellation of all the available CFD anchor points. 
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shown. 
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Mach 

Q, Bars 

Alpha, degrees 




1.30 

0.000001 

147 

4.00 

0.00001 

148 

6.00 

0.0001 

149 

10.00 

0.001 

150 

12.50 

0.010 

151 

15.00 

0.030 

152 

17.00 

0.065 

153 

19.00 

0.083 

154 

20.40 

0.100 

155 

21.40 

0.125 

156 

23.00 

0.150 

157 

25.00 

0.175 


26.00 

0.200 


27.00 

0.250 


29.00 

0.290 


30.00 

0.325 


31.00 

0.350 


32.20 

0.400 


34.00 

0.475 


36.00 



37.50 



40.00 



42.50 



45.00 



50.00 




Table. 2 Mach, Dynamic Pressure, and Angle of Attack conditions contained within the anchored CBAERO laminar and 
turbulent databases. The final tables are ‘square’ and thus contain 24 x 20 x 11 = 5280 total cases. 
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Fig. 23 “Two dimensional” tessellation of all the available CFD anchor points with the data base points (in red) shown. 



Fig. 24 “Two dimensional” tessellation of all the available CFD anchor points with the 5225 data base points (in red) 
shown in the full three dimensional Mach- Q- Alpha space. 
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Centerline Slice for M: 28.409700, Rho: 0.000350, T: 249.770004. AoA: 152.000000 Pressure vs. Z 

CFD Case: 29 

CF D File: . iCF D_Data/dpl r_h i_speed/cev. baseline .9kms.a28. lam 



Z, m 

Fig. 25 Mach 28.4, 152.0° Angle of Attack, Centerline Pressure Distribution - Anchored CBAERO solution compared with 
DPLR results. All of the available CFD points were used in the anchoring process. 


CenterlineS ce for M 28.409700 Rho 0.000350 T 249.770004 AoA 152.000000 Convective Qdot vs Z 

CFD Case: 29 

C F D File : . . /CF D_Data/d pi r_h i_speed/cev . base! i ne .9kms„a28 . la rn 



z__ 

Fig. 26 Mach 28.4, 152.0° Angle of Attack, Centerline Laminar Convective Heating Distribution - Anchored CBAERO 
solution compared with DPLR results. All of the available CFD points were used in the anchoring process. 
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Centerline Slice for M: 28.409700, Rho: 0.CCC35C. T 249.770004, AoA: 152.000000 Temperature vs. Z 

CFO Case: 29 

CFD File: . ./CFD_Datafdplr_hi _speed.xev.base ire.9kms.a28 Jam 



Z. nn 

Fig. 27 Mach 28.4, 152.0° Angle of Attack, Laminar Temperature Distribution - Anchored CBAERO solution compared 
with DPLR results. All of the available CFD points were used in the anchoring process. 



Fig. 28 Mach 28.4, 152.0° Angle of Attack, Centerline Laminar Shear Stress Magnitude - Anchored CBAERO solution 
compared with DPLR results. All of the available CFD points were used in the anchoring process. 
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Centerline Slice for M: 28.403700, Rho: 0.000350, T: 249.770004, AoA: 152.000000 Radiative Qdot vs. Z 

CFD Case: 29 

C FD File : . . /CFD_D ata/d plr_hi_speed/ce v. ba seline 9kms ,a28. 1 am 
CFD File: . . /C F D_Data/d plr_h i_s peed/neq airjresults/TC 1 5Lam.Rad 



z. m 

Fig. 29 Mach 28.4, 152.0° Angle of Attack, Centerline Radiative Heating - Anchored CBAERO solution compared with 
NEQAIR results. All of the available CFD points were used in the anchoring process. 



Fig. 30 Mach 28.4, 152.0° Angle of Attack, Surface Laminar Convective Heating - Anchored CBAERO solution compared 
with DPLR results. All of the available CFD points were used in the anchoring process. 
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Fig. 31 Mach 28.4, 152.0° Angle of Attack, Surface Radiative Heating - Anchored CBAERO solution compared with 
DPLR results. All of the available CFD points were used in the anchoring process. 
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